2207.02142v1 [astro-ph.GA] 5 Jul 2022 


arX1V 


DRAFT VERSION JULY 6, 2022 
Typeset using IATRX twocolumn style in AASTeX63 


Observation of $4716 - A star with a 4 year orbit around Sgr A* 


FLORIAN PEIGKER,’ ANDREAS ECKART,”’*? MICHAL ZAJACEK,? AND SILKE BRITZEN7 


17. Physikalisches Institut der Universitdt zu Kéln, Ziilpicher Str. 77, 50937 Kéln, Germany 
? Max-Plank-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
3 Department of Theoretical Physics and Astrophysics, Faculty of Science, Masaryk University, Kotldiské 2, 611 37 Brno, Czech Republic 


(Received January 1, 2018; Revised January 7, 2018; Accepted July 6, 2022) 


Submitted to ApJ 


ABSTRACT 


The ongoing monitoring of the Galactic center (GC) and Sgr A*, the central supermassive black hole, 
produces surprising and unexpected findings. This goes hand in hand with the technical evolution of 
ground- and space-based telescopes and instruments, but also with the progression of image filter 
techniques such as the Lucy Richardson algorithm. As we continue to trace the members of the S- 
cluster close to Sgr A* on their expected trajectory around the supermassive black hole, we present the 
finding of a new stellar source, which we call $4716. The newly found star orbits Sgr A* in about 4.0 
years and can be detected with NIRC2 (KECK), OSIRIS (KECK), SINFONI (VLT), NACO (VLT), 
and GRAVITY (VLTI). With a periapse distance of about 100 AU, $4716 shows an equivalent distance 
towards Sgr A* as $4711. These fast-moving stars undergo a similar dynamical evolution, since $4711- 
54716 share comparable orbital properties. We will furthermore draw a connection between the recent 
finding of a new faint star called $300 and the data presented here. Additionally, we observed a blend 
star event with $4716 and another new identified S-star $148 in 2017. 
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1. INTRODUCTION 


Almost 5 decades ago, Lynden-Bell & Rees (1971) 
proposed the existence of a supermassive black hole 
(SMBH) in the centre of our Galaxy. Following the 
claim by Lynden-Bell & Rees (1971), several generations 
of telescopes and instruments have revealed detailed in- 
formation about the structure and components of the 
center of our galaxy. For example, observations with 
NACO (Lenzen et al. 2003; Rousset et al. 2003) that 
were mounted at the Very Large Telescope (VLT) have 
provided access to precise stellar measurements of single 
stars on their relativistic (S2, see Parsa et al. 2017) and 
non-relativistic (SO-102, see Meyer et al. 2012) orbits 
around the central mass. This central mass can be de- 
scribed as the sum of the enclosed mass consisting of the 
compact mass (i.e. SMBH) and the extended mass (for 
example, nondetected faint stars, compact remnants, 
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and dark matter). 

From a historical point of view, Eckart & Genzel (1996) 
identified a few dozen stars orbiting Ser A* and named 
the related stellar members S-stars. These S-stars be- 
long to the S-cluster and can be characterized by veloc- 
ities up to several thousand km/s, with some of them 
moving on highly eccentric orbits (Ali et al. 2020). It is 
a remarkable achievement on both the technological and 
scientific level to analyze stellar orbits that pass around 
Ser A* at distances comparable to the size of our solar 
system. Observations underline the multifaceted (star- 
formation) history of the Nuclear Star Cluster (NSC) 
and the embedded S-cluster (Lu et al. 2013; Schédel 
et al. 2020; Nogueras-Lara et al. 2021). 

The authors of Eckart & Genzel (1996) analyzed the S- 
stars using the Lucy-Richardson algorithm (Lucy 1974) 
that stands for a high-pass filter. With these filtered 
data, Schddel et al. (2002) was able to analyze for the 
first time the orbit of $2, the brightest member of the 
S-cluster with an orbital period of about 16 years. Fur- 
thermore, Ghez et al. (1998) and Eckart et al. (2002) 
suggested that the compact and variable radio source 
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Sgr A* (Balick & Brown 1974) can be associated with an 
SMBH. The unexpected finding of young O/B stars close 
to Ser A* resulted in the formulation of the "Paradox 
of Youth" (Ghez et al. 2003) highlighting the unique- 
ness of the cluster. Due to the gravitational potential 
of Sgr A* but also the high dust temperatures of about 
200 K (Cotera et al. 1999) any possible star formation 
process is hampered (Jeans 1902). 

Due to a sensitive application of the high-pass filter in 
combination with a mild background subtraction, we are 
today able to observe stars that are even on shorter or- 
bits than S2 (magx ~14, see Rafelski et al. 2007) and 
S0-102/S55 (magx +17, see Meyer et al. 2012). These 
stars are S62 and $4711 with orbital periods of 9.9 years 
(magx 16, see Peifker et al. 2020a, 2021b) and 7.6 
years (magx #18, see Peifsker et al. 2020d), respec- 
tively. A natural by-product of these short-orbit stars 
is the amount of the minimized enclosed mass, which 
results in higher precision regarding the calculation of 
the compact and extended mass. Due to the confusion 
and blend stars (Sabha et al. 2012), the observation of 
S62 and $4711 at a distance of a few mas from Sgr A* 
can be highly challenging. On the other hand, stars at 
the edge of the S-cluster such as $29 are also not free 
of confusion. On their way towards the inner region of 
the S-cluster, the star interferes with several identified 
and non-identified S-stars causing a confused detection. 
From an observational point of view, the stellar con- 
fusion caused by close-by stars covers just one aspect. 
Due to the high extinction in the infrared of Ay =1-—3 
(Schédel et al. 2010; Fritz et al. 2011; Gautam et al. 
2019) towards the Galactic center, a photometric anal- 
ysis can be challenging, especially on large structures 
such as the NSC (Schédel, R. et al. 2009). However, 
substructures like the 40 mpc wide S-cluster should not 
be affected by extinction variations. The authors of Lu 
et al. (2013) for example use the common NIR K-band 
value Ax = 2.7 for their analysis of the entire S-cluster. 
However, the motivation for the observation of S-stars 
is surely related to the supermassive black hole Sgr A* 
and how stellar orbits evolve in the strong-gravity regime 
(Parsa et al. 2017; Hees et al. 2017, 2020). With a pre- 
cise knowledge of the orbital parameters, one is able to 
determine properties of Sgr A* which are necessary for 
theoretical models (see, e.g., Zajacek et al. 2020; Tur- 
sunov et al. 2020; Sukova et al. 2021). For example, the 
relativistic prograde pericenter advance of a star per its 
orbital period around a central mass can be determined 
by 


where M is the mass of Sgr A* (see Weinberg 1972). 
The gravitational length scale is defined by 


GM 
[= = 1.96 x 107-4 mpc ~ 4.9pas, (2) 


with M = 4.1 x 10°Mo, G the gravitational constant, 
and c the speed of light (see also Rubilar & Eckart 2001). 
The semi-major axis a and the eccentricity e are inferred 
from the orbital solution of a given star. Furthermore, 
the relativistic parameter is given by [ = = where 
rg represents the Schwarzschild radius of Sgr A* with 
rg = 2GM/c? = 1.21 x 10!7(M/4.1 x 10°Mo) cm and 
rp represents the pericenter distance of the stellar probe 
with 

Tp = a(l—e). (3) 


Comparison of relativistic parameters in units of 107+ 
for S2 (6.8), S62 (46), S4711 (5.6), and $4714 (64) 
(Peifker et al. 2020d) in relation to their related K- 
band magnitude reveals the challenging framework of 
the observations. Although S62 and $4714 are surely 
more suitable to investigate the relativistic influence of 
Ser A*, the faint emission in combination with the influ- 
ence of the dominating point spread function (PSF) of 
S2 hinders a confusion-free detection. Moreover, com- 
mon blend star scenarios as well as crowding problems 
caused by nearby/background stars are increasing the 
uncertainty range. One has furthermore take into ac- 
count a stellar density of over N = 10.0/arcsec? in a 
distance of 40 mpc (i.e., the size of the S-cluster) for 
sources with a K-band magnitude of mx = 17.5 — 18.5 
(Gallego-Cano et al. 2018). 

By inspecting the results of Gillessen et al. (2009b), 
Meyer et al. (2012), Parsa et al. (2017), Gillessen et al. 
(2017), and Peifker et al. (2020a,d) we come to the con- 
clusion that using the iterative Lucy-Richardson pro- 
cess is the most reliable way of analyzing data of a 
crowded area, such as the S-cluster. For example, the 
well-investigated orbit of $55 /S0-102 can be revealed us- 
ing a high-pass filter, such as the Lucy-Richardson al- 
gorithm, because overlapping PSF wings suppress star 
emission. 

This approach can be applied to most of the S-cluster 
members. Most recently, the orbit of the star S29 gained 
attention because of a disputed detection between 2015 
and 2021 (Peifker et al. 2020a; GRAVITY Collabora- 
tion et al. 2021; Peifker et al. 2021b). This confusion 
resulted in the observation of the periapse of a star in 
2021 that was interpreted as $29 by Gravity Collabo- 
ration et al. (2022a). Since we already showed the tra- 
jectory of S29 in Peifker et al. (2021b) on its trajectory 
towards Sgr A*, we will present an alternative inter- 
pretation with the introduction of a new S-star 54716. 
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We chose this name to follow up on our nomenclature 
presented in Peifker et al. (2020d) where we introduced 
several new stars, which we call $4711-S4715. These 
stars are orbiting Sgr A* on highly eccentric orbits with 
periods close to or less than 10 years. For $4716, we de- 
rive an orbital period of only 4 years. With a distance of 
about 12 mas during its periapse, $4716 has the shortest 
known orbital period around a supermassive black hole 
to date. 

The paper is structured as follows. In Sec. 2, we 
will provide a brief overview of the data used (see also 
Appendix A) and the methods. We furthermore com- 
pare different analyzing techniques. Then we present 
the results of our analysis in Sec. 3 and discuss the 
results in Sec. 4. Our final conclusions of the here pre- 
sented data can be found in Sec. 5. In the Appendix, 
we list the used data, provide further results, and pro- 
pose a guideline for the detection of single stars in a 
crowded stellar field. In addition , we will introduce 
another new S-cluster member called $148. We will ob- 
serve the emergence of a blend star event that justifies 
the checklist at the end of this work. In the following, 
we will focus on six stars: $148 (this work), 54716 (this 
work), 529 (see Peifsker et al. 2021b), S62 (Peifker et al. 
2020a), S29cGitlessen (See Gillessen et al. 2017), and S29 
(see Gravity Collaboration et al. 2022a). See the finding 
chart in Fig. 1 where we indicate the position of the lat- 
ter three in the white dashed circle. We will abbreviate 
Gillessen et al. (2017) by G17, GRAVITY Collaboration 
et al. (2021) with GColl2021a, and Gravity Collabora- 
tion et al. (2022a) with GColl2022a. Since we are aware 
of possible confusions for the reader, all abbreviations 
used in this work are listed in Table 1. 


2. DATA AND ANALYSIS 


Here we provide an overview of the data reduction 
and a short introduction of the used instruments. In 
addition, we list the tools that we have used to obtain 
our final images. We will furthermore describe the ap- 
plied background subtraction because of the dominating 
background in the S-cluster, especially close to Sgr A*. 
The used data is listed in the Appendix A. 


2.1. SINFONI 


The Spectrograph for INtegral Field Observations in 
the Near Infrared (SINFONI) was mounted up until 
2019 at the Very Large Telescope (VLT) at the European 
Southern Observatory (ESO) (Eisenhauer et al. 2003; 
Bonnet et al. 2004). It is now decommissioned. The 
instrument contained an Integrated Field Unit (IFU) 
which is responsible for the three-dimensional shape of 
the data (two spatial dimensions in x and y directions, 
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Figure 1. The K-band view of the Galactic center observed 
with NIRC2 (KECK) in 2019.30. This image is high-pass 
filtered and shows the position of several S-stars close to 
Sgr A* which indicated by a black x. The star in the white 
dashed circle shows a K-band magnitude of 16.3 mag while 
the star in the blue circle is 17.0 mag faint. Both K-band 
magnitudes are consistent with the reported values for S62 
by Peifker et al. (2020a) and S29 by Gillessen et al. (2017) 
and Peifsker et al. (2021b), respectively. Here the north is up 
and the east is to the left. 


and one spectral dimension along the z axis). With 
this setup, every pixel has a related spectrum. Con- 
tinuum images and line maps covering the H+K band 
can be extracted from the data cube. In this work, we 
limit the positional analysis to K-band continuum im- 
ages (2.0 — 2.2 um). 

For the data reduction, we used the ESO pipeline 
with DARK-, FLAT-, LINEARITY-, WAVE-, and 
DISTORTION-correction files. The final single 3d data 
cubes are then cropped to exclude non-linear edges and 
shifted in an array to a specific position (x;, y;). This 
position is determined by a manually created finding 
chart. Single data cubes (64 x 64 pixel) have a spa- 
tial pixel scale of 12.5 mas resulting in a field-of-view 
(FOV) of 0.8 x 0.8 arcsec with an integration time of 
600 seconds. The final data cubes have a FOV of about 
1.2 x 1.2 arcsec. The stacked integration time depends 
on the number of single data cubes. Adaptive optics 
(AO) for every observation is enabled and, if possible, 
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Star related Publication abbreviated Publication 
$148 this work - 
54716 this work - 
$29 Peifker et al. (2021b) - 
S62 Peifker et al. (2020a) = 
S29Gillessen Gillessen et al. (2017) G17 
$29 Gravity Collaboration et al. (2022a) GColl2022a 
GRAVITY Collaboration et al. (2021) GColl2021a 


Table 1. Abbreviations used in this work. We list the investigated stars in this work and the related publications that should 


be consulted by the interested reader. 


used with an optical natural guide star in the proximity 
of Sgr A*. 


2.2. NACO 


The Nasmyth Adaptive Optics System (NAOS) Near- 
Infrared Imager and Spectrograph (CONICA), abbrevi- 
ated to NACO, was decommissioned in 2019 at the VLT 
(Lenzen et al. 2003; Rousset et al. 2003). The imager 
operated in the H-, K-, L-, and M-band with an infrared 
wavefront sensor for AO. For the nonarchival observa- 
tions, we used a 5.0 x 5.0 arcsec box with a random 
dither pattern. For every here analyzed observation, the 
bright pulsating supergiant IRS7 located ~ 5.5” north of 
Ser A* is used for the AO correction. The final images 
are shifted and added to increase the on-source integra- 
tion time. The spatial pixel scale of the here analyzed 
data is 13.3 mas with a FOV of 1024 x 1024 pixels. 
The observations were executed in the K,-band that is 
covering 1.97 — 2.32um. 


2.3. NIRC2 


At the KECKII telescope (located at Hawaii/USA), 
the Near-Infrared Cam2 (NIRC2) is mounted (Matthews 
& Soifer 1994; Nelson 1997). The imager operates in the 
H-,K-, and L-bands using AO to correct for atmospheric 
turbulence. The FOV of single images is 1024 x 1024 
pixel with a spatial pixel scale of 9.9 mas. From the 
KECK Observatory Archive (KOA'), we downloaded 
the public data that was observed in the K,-band cov- 
ering the wavelength range 1.94 — 2.29um. 


2.4. OSIRIS 
Like SINFONI, the OH-Suppressing Infrared Imaging 
Spectrograph (OSIRIS) produces 3d data cubes with 2 
spatial and 1 spectral dimension (Larkin et al. 2006; 
Mieda et al. 2014). The KECKII mounted instrument 
uses an Integrated Field Spectrograph (IFS) with AO. 
The instrument itself also contains a science camera 
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that produces continuum images of the related scien- 
tific object. In this work, we will focus on the results 
of this science camera. The OSIRIS imager has a FOV 
of 2048 x 2048 pixel. Here, the presented data is ob- 
served with the Kn3 filter covering the K-band from 
about 2.12 um to 2.22 um with a spatial pixel scale of 10 
mas. We shifted and added the downloaded data from 
the KOA to increase the on-source integration time to 
about 8 seconds. 


2.5. High-pass filter and methods 


In general, high-pass filters are used to suppress back- 
ground noise and dampen the wings/airy rings of a PSF. 
Since objects can move through intersecting regions that 
are polluted by overlapping airy rings or wings, one can 
reduce this influence by a mild background subtraction. 
This process is not only needed for the analysis of a 
crowded region, but it maximizes the validation of the 
data. Since artefacts are removed by the background 
subtraction, only emission above the noise level is ob- 
served. Vice versa, if the image is dominated by noise, 
every stellar emission is suppressed and hence, no ob- 
ject can be observed. In this work, we use two dif- 
ferent techniques, namely the smooth subtraction (SM) 
and the Lucy-Richardson (LR) algorithm, with the fo- 
cus on the latter. While the SM-algorithm is a reliable 
tool, the resulting images still suffer from noise and ex- 
tended diffraction minima that can be observed espe- 
cially around $2. On the other hand, the LR-algorithm 
is sensitive to the chosen input parameters and should be 
complementary to the SM-algorithm to avoid a confused 
analysis. Even though we described the two filter tech- 
niques numerous times (see, e.g., Peifsker et al. 2020a,d, 
2021b), we will list the necessary steps for transparency 
in the following. In Fig. 2, we display K-band im- 
ages observed with NIRC2 (KECK) and NACO (VLT). 
We compare these continuum images with two differ- 
ent high-pass filtering techniques, namely, the SM- and 
LR-alogrithms. Figure 2 also shows the difficulties of ob- 
serving and finding new and unknown stellar S-cluster 
members. Hence, we will provide a checklist in Ap- 
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NACO, No filter 


NIRCAM2, No filter 


Smooth subtract Lucy Richardson 


Figure 2. Comparison of the K-band continuum emission with two different instruments and two high-pass filter techniques 
in 2010. The left NIRCAM2 image is adapted from Meyer et al. (2012) while all the other plots show NACO data. In the far 
right image, we indicate the orbit of S62 and in this work discussed new S-cluster member $4716 with lime- and blue-colored 
trajectories, respectively. Please compare the observation of $55 in every presented panel. In every image, north is up, east is 


to the left. 


pendix J that establishes a frame for the identification 
of unknown sources. 


2.5.1. Smooth-subtract algorithm 


As an advantage, the SM-filter is easy to use and han- 
dle. The final image I, is a convolution of point sources 
and a PSF. The final image I,;, however, is a combi- 
nation of all observed frequencies whereas a smoothed 
version of I) is representing a low-pass filtered image 
Tjow: To derive the final SM-filtered image, we can write 


Lan = dise = Tish (4) 


where Ipign is the desired result. With the above equa- 
tion, it becomes eminent that the only input param- 
eter of this procedure is the smoothing PSF that we 
denote with PSF smooth. As a general procedure, one 
set PSF smooth = PSFcauss where PSF causs represents 
a Gaussian filter with a size that matches the PSF of 
the data. As it is mentioned before, the subtraction of 
a PSF often results in a negative flux ring around the 
brightest source, making it impossible to analyze the 
close-by regions. 


2.5.2. Lucy-Richardson algorithm 


As a major keystone, the PSF and background sub- 
traction for the LR-filter determines the output and 
should be used as a complementary process or be cross- 
checked with different datasets. Nevertheless, there are 
different procedures available to extract a PSF (see for 
example, StarFinder, Diolaiti et al. 2000). We will fo- 
cus here on a manual approach to have full control over 


the parameters. As a first step, we check isolated stars 
in the crowded region. Here, a star is labeled isolated 
if the distance to the next near-by source is at least 
one PSF. Possible reference stars are, for example, 54, 
S6, S7, S30, S65 at a distance of 0.2-0.8 arcsec from 
Ser A* (Ali et al. 2020). With the selection of these local 
stars, we minimize anisoplanatic effects (Schddel 2010). 
Since this approach cannot always be fulfilled (for ex- 
ample, the SINFONI FOV is highly limited in size), at 
least S2 can be used to get the first idea of the true 
PSF. The parameters are derived by fitting different- 
sized Gaussians to the stars. Common parameters are 
the « (FWHM1) and y (FWHM2) values as well as the 
angle of the PSF. Another obvious parameter is the mag- 
nitude of the reference stars. With magx = 14.1 mag, 
S2 is the brightest S-cluster member and thus is the 
most prominent star in the data investigated. Hence, 
the reference stars should not be too faint, i.e., too close 
to the background since this could affect the shape of 
the PSF. A K-band magnitude of < 15.5 mag is a reli- 
able choice since the background emission accounts for 
about 2-5% of the peak flux of the reference source. The 
next step is the application of the LR-algorithm to get 
rid of the dominant side lobes of the image PSF. Af- 
ter this deconvolution step, the resulting delta function 
map is convolved with a PSF with the size dimensions of 
= (FWHM1 + FWHM2)/2. In the resulting convolved 
image, the brightest and isolated sources are again fit- 
ted with a Gaussian of different sizes if the image PSF 
is not circularly shaped. This process will be repeated 
until the fitted PSF converges to a circular shape with 
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Figure 3. Flowchart of iterative process for using the LR 
algorithm. The input image is a FITS file. However, other 
processable formats are allowed. The PSF selection is usually 
done outside the S-cluster because of the number of isolated 
sources (see text for details). A FOV can suffer from a vari- 
able PSF which increases the challenge of finding a suitable 
solution. However, we limit the analysis to a small FOV 
where we assume, that the PSF is stable to first order (i.e., 
variations of less than 0.25 pixel are allowed). Because S2 is 
usually not an isolated source, it should only be used if no 
other stars can be found (e.g., the SINFONI data). After ap- 
plying the LR algorithm with iterations of 10k, we convolve 
the resulting delta map with a GAUSSIAN. The convolved 
image is the solution if the image does not suffer from arte- 
facts. The process can be cross-checked by the position of, 
for example, 82. 


a minimized elongation. The methodology of this pro- 
cess is shown in Fig. 3. We note that the process de- 
scribed with the flow diagram shown in 3 only serves 
as an overview. If, for example, the overall quality of 
the data are reduced, this iterative will not result in the 
finding of a suitable PSF. 

Although background subtraction can help to find the 
real PSF of the image, a much more applicable approach 
is the construction of an Airy-ring filter that we call 


PSF,. Once the iterative process described here is com- 
plete, a suitable PSF is constructed which we denote by 
PSF pst. We create a new filter with 


PSF, ® PSF psp = PSFapst (5) 


where PSF asp denotes the convolved result of the input. 
Using the LR-algorithm with PSF asp¢ in combination 
with a background subtraction, we derive robust results 
showing stars at their expected orbital positions. We 
note that the process described is the same for every 
data set. 


2.5.3. Sophisticated background subtraction 


By analyzing the S-cluster, one is permanently con- 

fronted with the dominant and variable background 
(Tep et al. 2021). One example that underlines the 
complexity of the topic is the faint star $300 that is 
observed with GRAVITY (GColl2022a). The existence 
of $300 with its apparent K-band magnitude of almost 
20 mag adds to the cusp of faint stars that are under 
the detection limit regarding the imaging data observed 
with, for example, SINFONI, NACO, or NIRC2. 
As shown by Sabha et al. (2012), the diffuse background 
emission increases towards Sgr A*. Because of the struc- 
ture of the S-cluster, the increasing/decreasing back- 
ground is not symmetrical with respect to Sgr A*. For 
the moment, we will ignore the nonuniform shape of 
the cluster. Assuming an image without errors (e.g., 
dead/hot pixels or cosmic rays), we can define every 
pixel as 


signalops = signalop; 7 signalpackground (6) 


where signal,,; is related to an arbitrary object (e.g., a 
star). Naturally, signal,,; can be equal to zero and the 
related pixel information is limited to the background 
emission. Since the variable background is a function of 
the position r, it is 


f(r) = signalpackeround (7) 


and 
—grad(fi)) # 0 (8) 


where grad(f;,)) is directed towards Sgr A*. We define 
z= fir) = x-y and set R.A. = x, DEC = y. We use 


Of Of of 
_ A) = ( (2) yey = —(r)e, (9 
grad( i) (an (x), Or, (r)) d Or; (r)é ( ) 
and since we observe the on-sky projected view, we 
limit the above equation to 7 = z. The modulus 
yields y = i, ignoring negative values. Using now 


signal,pj ~ 0, we can investigate different scenarios for 
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Figure 4. A trend of the background inside the S-cluster. 
Here we normalize the size of the S-cluster to unity. Please 
note that this trend can differ depending on the direction 
because of the non-uniform shape of the S-cluster. Obser- 
vations of this trend are shown in Sabha et al. (2012) and 
Peifgker et al. (2021c). 


the background subtraction procedure. In Fig. 4, we 
show the trend for the diffuse background in an arbi- 
trary direction where the origin is defined as the posi- 
tion of Sgr A*. As we pointed out before, the shape of 
the S-cluster is nonuniform and hence, the trend in Fig. 
4 might differ for different directions. As a consequence, 
we introduce a variable parameter “b” with b € Rt/{0} 
and y = 4 -b. We set b = 1 for the example in Fig. 4 
and randomly disperse 5 stars with different magnitudes. 
These represent different background subtraction cases: 


e star 1: The star is above the background and 
can be observed. If, however, star 1 moves closer 
towards Sgr A*, the stellar magnitude might in- 
crease. For a confusion-free detection, a back- 
ground subtraction can be applied. One should 
note that an adaptive background subtraction 
could result in a decreased chance to observe star 
5. The stellar object star 4 can be observed re- 
gardless of the background subtraction. 


star 2: This star cannot be observed at this dis- 
tance but contributes to the variable background 
emission. It can be assumed that this star also 
orbits Sgr A* and creates fluctuations in the back- 
ground or a blend-star scenario. It is important to 
note that this star can be observed if the defined 
orbit allows a sufficient distance with respect to 
Ser A*. 


e star 3: The star is confused with the background 
and cannot be observed at this distance. It can be 
treated like star 2. 


star 4: With a magnitude above the highest noise 
level, this star can be observed regardless of the 
background or the distance to Sgr A*. Of course, 
if another star with a comparable magnitude is 
close by, the confusion is increased. 


star 5: Observing star 5 is possible but challeng- 
ing because of the prominent noise (e.g., $55 and 
S62). Again, an increased distance with respect to 
Ser A* enhances the probability to observe a star 
like star 5. 


We refer the interested reader to the Appendix of 
Peifsker et al. (2020a) where we consider additional as- 
pects of the background emission and source crowding. 


2.6. Methodology 


Here, we outline the process of extracting data points 
from our data set. Since we use data from 5 different 
instruments that are mounted at VLT and KECK, we 
have to take into account possible small offsets between 
the different reference frames. The general problem 
is described Gillessen et al. (2009a) where the authors 
show small offsets between the derived orbital solution 
for $2/S0-2 using KECK and VLT data. In contrast, 
we cannot confirm the given offsets in Gillessen et al. 
(2009a) of Ax = —3.7+0.6 mas and Ay = —4.1+0.6 
mas. Using a 17 year data baseline to derive orbital so- 
lutions for S-stars like, for example, S29 observed with 
NACO and SINFONI, we find no significant offset. As 
we show in this work, the orbital solutions are inter- 
changeable. As discussed in Peifker et al. (2021c), the 
background noise is the most dominant contributing fac- 
tor to uncertainties. 

Matching Keplerian stellar solutions covers one aspect 
of the data treatment. In addition to noise (caused by 
the background or the instrument), the exact position 
of Sgr A* is the another strong source of uncertainties. 
As described in Gillessen et al. (2009b) and Parsa et al. 
(2017), the SiO maser stars in combination with a linear 
transformation is a suitable way to derive the position 
of Sgr A*. This is reflected in a precise determination of 
the orbital elements of S52 allowing to investigate, for ex- 
ample, the Schwarzschild precision. In agreement with 
Gravity Collaboration et al. (2018), the authors of Parsa 
et al. (2017) demonstrated that the NACO data provide 
a high level of precision. Furthermore, the general pro- 
cess for deriving the position of Sgr A* in the GRAVITY 
and NACO data is comparable and results in the match- 
ing outcomes as the authors mentioned above find a rel- 
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ativistic parameter of 0.88 (Parsa et al. 2017) and 0.9 
(Gravity Collaboration et al. 2018). The VLT results are 
accompanied by the analysis of Do et al. (2019), where 
the authors estimate a relativistic parameter of 0.88 for 
S0-2/S2 observed with the KECK telescope. Therefore, 
we can conclude that the well-observed orbit of S2 al- 
ready includes several detector and reference frame cor- 
rections. Hence, we use the orbital elements of S52 from 
Do et al. (2019) and GColl2022a for the analysis of the 
KECK and VLT data, respectively. 


3. RESULTS 


Here we present the results of our analysis. We in- 
troduce $4716, a new member of the fast and faint 
S-cluster population (FSS-population, see also Peifsker 
et al. 2020d). 


3.1. The Keplerian orbit of $4716 


As we constantly analyze the existing and new imag- 
ing data of the GC, we inspected the data presented 
in Peifker et al. (2020a), Peifker et al. (2020d), and 
Peifsker et al. (2021b) with the here presented fine-tuned 
approach. Based on this analysis, we recognized an 
unidentified source between Sgr A* and $27. Especially 
the data of 2007 and 2019 showed a bright source that 
does not fit in our stellar orbit list (Ali et al. 2020). Be- 
tween 2003 and 2020, we identify this stellar candidate 
in 16 epochs (see Table 8, Appendix A and Fig. 5) for a 
selected overview). Limitations are the stellar interfer- 
ence from various sources such as, for example, $2, $55, 
and $62. This can be observed, for example, in 2004, 
2011, and 2012 where S62 (Peifker et al. 2021b) com- 
plicates a confusion-free observation®. In addition, the 
configuration of the Fine-tuned Artificial PSF (FAPSF) 
and the intensity of the related wings (i.e., Airy rings) 
are sensitive input parameters to the Lucy Richardson 
algorithm. As mentioned before, a list to verify the find- 
ing of FSS-stars is presented in Appendix J. 

Using the positions from the NACO, SINFONI, NIRC2, 
and OSIRIS data analysis results in the Keplerian so- 
lution for the stellar candidate (named $4716) with an 
orbital period of 4 years (please find the R.A.-time and 
Dec.-time plot in Fig. 6). Recently, GColl2022a ob- 
served in detail the periapse of $29 that coincides with 
the here presented Keplerian solution for 54716. We 
will comment on this confusion in the Discussion sec- 
tion. However, we include the positional GRAVITY 


2 Please consult Peifker et al. (2020d) and more recently, Peifker 


et al. (2021b). 


3 See the subplot for 2010.2 in Fig. 5. The star on the right next 


to $4716 is S62. 


data of 2021 for S29 as presented in GColl2022a and 
found a remarkably good match with the orbit of $4716 
that is based on SINFONI and NACO observations cov- 
ering the epochs 2003-2020. In Table 2, we present the 
orbital parameters of $4716. We additionally list the 
orbital elements of S2 since this well-observed star is 
used to determine the position of Sgr A* in agreement 
with the results described in Parsa et al. (2017). The 
authors of Parsa et al. (2017) linked the infrared posi- 
tional frame on a few milliarcsecond scale to the posi- 
tions of known GC MASER stars obtained from the Very 
Large Array (VLA) radio observations. The orbital ele- 
ments listed for $2 in Table 2 are taken from Do et al. 
(2019) and agree reasonably well with the VLA results 
mentioned but also with GColl2021a. We find a reason- 
able agreement of the estimated Keplerian orbit with the 
maximum likelihood given by the Markov Chain Monte 
Carlo (MCMC) statistics (Table 2). While we limited 
the range of the pericenter passage to 2017.42 + 0.2 
years, we allowed every possibility for the mass and dis- 
tance of Sgr A*. We will elaborate on the statistical 
results for the mass and distance of Sgr A* in Appendix 


C. 


3.2. Parameters of $4716 


The analyzed epochs in this work are showing $4716 
as an isolated object. In these epochs, the identifica- 
tion resists the influence of the dominant PSF of 82. 
However, the variable footprint and distance of the $2 
PSF affect the determination of the brightness of $4716. 
Therefore, we averaged the K-band magnitudes for indi- 
vidual identifications between 2003 and 2020 (see Table 
8, Appendix A) to cover a large data baseline and coun- 
terbalance the mentioned fluctuations. Analyzing the 
data shown in Fig. 5, we find an average K-band mag- 
nitude of 17.02 + 0.22 mag (Fig. 7) by using 


magsa7ie = —14.15 + 2.5 x log(countsg4716 /countsg2) 

(10) 
where S2 is used as the reference source (Schédel et al. 
2002). Furthermore, we use 


fsa7ie = fs2 X 190~0-4(magsa7i6 — mags2) (11) 


which is adapted from Sabha et al. (2012) to derive a K- 
band flux of fg471g¢ = 1.0340.15 mJy with fso = 14.72 
mJy and magsgg = 14.1 (Schédel et al. 2002). Using the 
estimated magnitude, we follow the analysis of Peifsker 
et al. (2020d) to calculate the mass of $4716 with 


Msa716 
(Mss 


log = k-magga7ig + 5 (12) 


© 


where k equals -0.1925 and b = 3.885 (see Appendix 
I). With Eq. (12), we find a stellar mass for $4716 of 
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Figure 5. Timeline of the new identified star $4716 on its Keplerian orbit around Sgr A*. Here we show 12 epochs between 
2003 and 2019 in the direct vicinity of Sgr A*. The data of 2006.3, 2019.3, and 2019.5 are observed with NIRC2/KECK. In 
2013.4, 2014.5, and 2015.3, we show observations executed with SINFONI/VLT. All other epochs are observed with NACO/VLT. 
Ser A* is indicated by a cross with $4716 at the expected position marked with a green circle. We want to underline that for 
every shown epoch, S62 can be found at the expected position based on the orbital elements presented in Peifsker et al. (2020a) 
(except for 2015 where S62 is confused with S2). Therefore, we present exactly the same figure in the Appendix with the orbit 
of S62. Here, north is up and east is to the left. Every panel shows a 20x20 pixel image of the direct vicinity of Sgr A* and 
is normalized to the peak emission of S2 (which corresponds to a K-band magnitude of 14.1). The contrast is choosen to show 
also fainter sources that may seem brighter than the estimated magnitude. 


Mgarie = 4.04¢7 Mo. Using the relation L ~ M&y74¢ 
with a = 3.5, we estimate a stellar luminosity of 
logts41s = 2.127940. 

With the semimajor axis of $4716 listed in Table 
2, we find with Eq. (1) a relativistic Schwarzschild 
precession of about 0.25° ~ 15.1’ and a periapse 
distance of r, = 0.47 — 0.48 mpc & 97 — 99 AU. 
Furthermore, we estimate a relativistic parameter of 
[ =rg/rp ~ 8.0—8.2 x 1074 which is comparable to 
that of S2 (6.8 x 10~*, see Peifker et al. 2020d). 


Due to the isolated position of $4716 in 2007, we ex- 
amine the related spectrum (Fig. 8) and find a Doppler- 
shifted Bry emission line at 2.1625 um. Using the Bry 
rest wavelength at 2.1661, we derive an LOS velocity of 
about vgry = —430+55 km/s. In addition to the blue- 
shifted Bry line, we detect two Hel absorption lines. The 
helium line with the transition 2p'P° — 2s!S at a rest 
wavelength at 2.058m shows the same blue-shift as the 
Bry line with an LOS velocity of vger = —437 + 55 
km/s. The second detected Hel absorption line forms a 
shoulder feature (Habibi et al. 2017; Peifker et al. 2020d) 
that can be observed for many O/B-stars in the sur- 
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Source a (mpc) e i(°) w(°) (°) ldisceat( VE) 
52 5.09 + 0.01 0.886 + 0.002 133.49 + 0.40 65.89 + 0.75 227.46 + 1.03 2018.37 + 0.02 
S4716 (estimated) 1.93 +0.02 0.7562 0.02 161.2424 2.80 0.073 40.02 151.5441.54 2017.41 + 0.004 
S4716 (max. likelihood) 1.94+ 0.02 0.756+0.02 161.13 + 2.80 2.25 40.02 153.55 41.54 2017.41 + 0.004 


Table 2. Orbital elements of $2 and $4716. The estimated values for $4716 are based on a Keplerian fit model while the 
maximum likelihood stems from the MCMC simulation. See Appendix B for a discussion of these values. 
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Figure 6. Orbit of $4716. Here we show the R.A. and Dec. 
component as a function of time. The epoch ¢ is indicated 
by the numbers next to the data points. Because of the 
limited space, we choose the format "¢- 2000" for the nu- 
merical values. As indicated in this figure, we identify $4716 
in 16 epochs. In additon, we incorporate the observations 
by GColl2022a (see text for details). Blending events and 
the superposition with S2 can not be excluded and should 
be taken into account. For a detailed list of the data points, 
we refer to Table 8, Appendix A. 


rounding S-cluster. Since the spectrum of $4716 (Fig. 
8 and Fig. 18, Appendix D) is comparable to $4711 
(Peifker et al. 2020d), we assume the same spectral type 
for both stars. Therefore, using the atlas of Hanson et al. 
(1996) and the approach of Peifker et al. (2020d), we 
find similarities with a B8/9-V star. We follow up on 
this result in the following sections. This LOS velocity 
is consistent with the orbit fit result presented in Fig. 
6. The final LOS fit is shown in Fig. 17, Appendix D. 
Compared with the periapse LOS velocity of about 3000 
km/s, the star is slow during its apoapsis, indicating the 
possibility of detecting $4711. The projected on-sky plot 
of the orbit and the implemented data points displayed 
in Fig. 9 underline the probability of the observation of 
54716 during apoapsis. 

Due to the high blue shift and width of the Bry- 
line, the detection of the Hel doublet 77F — 4°D and 
7'F — 4'D is confused. Since the shape and velocity 
of the Hel lines at 2.161 um and 2.162 um are also in- 
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Figure 7. K-band magnitude of 54716 between 2003 and 
2020. The red dashed line represents the first order poly- 
nomial fit. The blue points show the data. We derive an 
average K-band magnitude of 17.02 + 0.22 mag. The mag- 
nitude is well inside the lo uncertainty range. 
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Figure 8. Spectrum of $4716 extracted from the SINFONI 
data cube of 2007. The rest wavelength of Hel 2p'P° — 2s'S 
at 2.058um is indicated with a dashed line. At 2.161yum and 
2.162um, the Hel doublet is located. The prominent absorp- 
tion feature of the Bry line at 2.1624m is blue-shifted with 
respect to the rest wavelength at 2.166um. The redshifted 
spectrum of $4716 observed in 2009 is shown in Fig. 18, 
Appendix D. 
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Figure 9. Projected on-sky view of the orbit of $4716. The 
red circled data points represent the GRAVITY measure- 
ments in 2021 are taken from Gravity Collaboration et al. 
(2022a). Please see Table 8, Appendix A for an overview of 
the shown data points. 


fluenced by the rotational speed vsin(i) and magnetic 
activity, we can only state a lower limit of the Doppler- 
shift of Vuel,,, = —210 km/s. Because this estimate is 
influenced by the rotation of the star (and hence not nec- 
essarily completely correlated with upoppier), the lower 
value is in reasonable agreement with the former de- 
rived blue-shifted velocities of about 430 km/s. Fur- 
thermore, the difference between the 7!F — 4!D line 
with the rest wavelength at 2.162 4m and the spectral 
position of v#er,,,, 18 1.5 x 10-34m. The difference is 
equal to 15. Now, we assume that the FWHM equals 
two times the Dopplershift difference between the line 
emission and the rest wavelength. This assumption is 
not entirely true due to the overall structure of the con- 
fused Hel doublet lines (Fig. 8). Therefore, we use the 
derived lower limit of vHer,,,, = —210 km/s and the up- 
per limit for the width of the line 2 x vyer,,, = —420 
km/s to give a rough estimate of the stellar rotational 
velocity Vrot = vsin(i) where we set v = 315 km/s 
based on previous assumptions. However, this results in 
Urot = 288+ 96 km/s. 


3.3. Stellar type 


Here, we follow up on the stellar properties derived 
in the former subsection. We use the luminosity of 
log “4z18 21 and the mass with Ms47ig = 
4.04+7 Me of $4716 as initial conditions for the BONN- 
SAI (BONN Stellar Astrophysics Interface)* (Schneider 


4 The BONNSAI web-service is available at www.astro.uni-bonn. 


de/stars/bonnsai. 


et al. 2014) analysis. The free available web-service 
BONNSAI is used to compare observed parameters with 
stellar models (see Brott et al. 2011). As a boundary 
condition, the star analyzed should have a related mass 
between 3-100 Mo. The motivation of using the ap- 
proach is the access to parameters that are not observed 
or to verify parameters, that suffer from increased confu- 
sion due to, for example, noise. However,we exclude v,o4 
as a free parameter to investigate if BONNSAT finds a 
value that is comparable to the one estimated in the for- 
mer subsection Sec. 3.2. Although our number of priors 
is low, the resulting outcome is in good agreement with 
the input values (Table 3). We find an effective temper- 
ature of about T.¢¢ = 12500 K which underlines the 
classification of a B8/9-V star. Despite the wide un- 
certainty of the stellar rotational velocity, the outcome 
agrees with 1 o with our independently derived value in 
the previous section. The result of the stellar rotational 
velocity from the BONNSAI simulation marks an inde- 
pendent control point and confirms our assumptions in 
Sec. 3.2. 


3.4. Mass and distance of Sgr A* 


Since $4716 with its short and small orbit offers a suit- 

able choice of pinpointing to an exact value, we derive 
with the help of MCMC simulations the most likely dis- 
tance and enclosed mass of Sgr A*. For each plot, we 
use 10k iteration steps and open boundaries for the mass 
and distance ensuring a nonbiased outcome. 
The presented results in Fig. 14-16 (see Appendix C) are 
used to estimate an averaged mass of (4.023 + 0.087) x 
10° Me for Sgr A* where the uncertainty represents 
the mean absolute error. Considering the distance of 
Ser A*, several attempts have been made to provide a 
consistent value. With the recent publications of Do 
et al. (2019) and Gravity Collaboration et al. (2022b), a 
4% discrepancy is introduced regarding the distance of 
Ser A*. Here we incorporate both results because the 
averaged distance of the MCMC statistics for the $4716 
orbit gives (8.028 + 0.199) kpc or alternatively, 8.028 
kpc + 2%. 


3.5. Source identification in 2020 


As it is shown in Ali et al. (2020), the number of S- 
stars with a related Keplerian orbital solution is much 
higher as discussed in Gillessen et al. (2009b) and G17. 
We have shown in several publications, that the confu- 
sion in the S-cluster is enhanced because of nonidenti- 
fied stellar objects (see $148 in this work and Fig.2 in 
Peifsker et al. 2021b). With this degree of confusion, 
the observation of S62 or $4711 becomes challenging. 
However, gap-free long-term surveys of the vicinity of 
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Table 3. Results of the BONNSAI simulation of stellar evolution. 


Ser A* are necessary to archive a uniform overview of 
the S-cluster. Since NACO and SINFONI were decom- 
missioned in 2019, we used KECK data to continue the 
analysis of the S-stars with data that provides a com- 
plete overview of the cluster. In Fig. 10, we show the 
data that was observed with the OSIRIS science cam 
with a spatial pixel scale of 10 mas. Consistent with 
the analysis covering 2003-2019, we find $4716 at the 
expected position in 2020.42 and 2020.58 (see also Fig. 
6 for the related orbit position). As we follow the mo- 
tion of S-stars (Ali et al. 2020), we do find the listed 
traced stars at the expected positions. With the orbital 
elements for S29@inessen/529 given in Gillessen et al. 
(2009b), G17, and GColl2021a, we do not find a stel- 
lar source matching the high-precision measurements in 
the related publications. 

However, it seems that S42 comes close to the orbital po- 
sition predicated by the authors mentioned above. Due 
to the increased crowdedness in the vicinity of Sgr A*, 
this star could be associated with $29 and would fit the 
orbital position derived by G17 using a 25-year data 
baseline. However, the derived pericenter passage for 
$29 is in 2025, excluding the possibility that GColl2022a 
observed this star during its closest approach in 2021. 
Additionally, we note the observation of $4711 as part of 
the FSS population in 2020 (Peifker et al. 2020d). This 
stellar S-cluster member can be found at the expected 
and predicted orbital position (see Fig. 10). 


3.6. Orbital properties of $4716 and spin of Sgr A* 


The orbit of $4716 is the most compact orbit within 
the S cluster. With the semi-major axis of only a = 
1.95 + 0.02 mpc = 402 +4 AU (given by the maximum 
likelihood value), its orbital period is Porp = 4.02 + 0.08 
years. The pericenter distance is rp = a(l —e) = 
(99+8) AU = (2494+ 204) r, and the apocenter distance 
isr, = a(1+e) = (706411) AU = (177834273) rg, which 
makes it the first stellar orbit around the SMBH that is 
fully inside 1000 AU, i.e. on the Solar-system length- 
scale. Given the eccentricity of e = 0.754 + 0.02, the 
velocity increases significantly by a factor of (1+e)/(1— 
e) ~ 7.13 from the apocenter value of (1116+53)kms7! 
to the pericenter value of (7956 + 381) kms~?. In other 
words, when scaling to the light speed, the orbital ve- 
locity of S4716 increases from (0.37+0.02)% of the light 
speed to (2.65+0.13)% of the light speed, which is within 


Table 4. Summary of relativistic parameters of 54716. 


Parameter Value 
Relativistic parameter I (8.02 + 0.68) x 10~* 
Schwarzschild precession 6¢ [arcmin] 14.84 1.3 
max. gravitational redshift [km s~*] 120 + 10 
max. transverse Doppler shift [kms~*] 106 + 10 
max. combined redshift [kms~*] 226 + 20 
LT precession rate [x = 0.5, arcsec yr *] Lie OS 


lo uncertainty comparable to the $2 star at the pericen- 
ter. 

Using the definition of the relativistic parameter 
(Parsa et al. 2017; Peifsker et al. 2020d), 


ec (13) 
Yp 
we calculate [ = (8.02 + 0.68) x 10~* for $4716, whose 
mean value is larger by a factor of ~ 1.2 and ~ 1.4 than 
for $2 and $4711 stars, respectively. 

The relativistic quantities of interest — the prograde 
apsidal (Schwarzschild) precession 6¢, gravitational red- 
shift vgr, and the Lense-Thirring precession Qur — are 
essentially proportional to the first power of T (44, ver) 
or its second power (Qur). Here we estimate these val- 
ues as they could be of interest for future monitoring, 
e.g., with the Extremely Large Telescope. 

The prograde relativistic precession (Schwarzschild 
precesion) may be estimated as 


61GM Bt 
ca(1—e?) 1+e 


which for $4716 can be evaluated as 6¢ = (14.8 + 1.3) 
arcmin. The positional shift of the apocenter can be es- 
timated as As = r,Ad¢d = 3raIl = 0.015 mpc ~ 0.4mas 
that corresponds to less than 5% of a pixel. 

Another general relativistic effect that affects the 
wavelength of the observed spectral lines is the gravita- 
tional redshift that for the pericenter may be evaluated 
as 


5b = 


iM (14) 


where the approximation on the right is valid for the 
pericenter. At the pericenter of $4716, the gravitational 
redshift amounts to vg. = (120+ 10)kms~*. The gravi- 
tational redshift cannot be practically disentangled from 
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Figure 10. Observation of $4716 in 2020.42 and 2020.58 with OSIRIS. We highlight the detection of $4711 at the expected 
orbital position (Peifsker et al. 2020d). On these scales and sensitivity, a star like $300 are not visible (GColl2022a). However, 
S300 is part of the variable background and serves as an example for the challenging analysis of the imaging data. See the text 


for details. North is up, east is to the left. 


the special-relativistic transverse Doppler shift, which 
for $4716 amounts to v%, = c[(1 — (vp/e)?)7/? — 1] = 
(106+10)kms7! at the pericenter. The combined max- 
imum relativistic redshift at the pericenter for $4716 
then is (226 + 20)kms7}. 

Finally, the compact orbit of $4716 is affected by the 
frame-dragging or the Lense-Thirring (LT) effect in the 
spacetime of a rotating (Kerr) black hole. The orbit- 
averaged change of the angular momentum direction can 
be calculated with 

dL 2G? M2 
(=) - al — €2)3/2q3 * xe (16) 
where y is the direction of the Sgr A* spin and L is the 
angular momentum direction of $4716. The frequency 
of the angular momentum precession can be evaluated 
as 
_ 2G? Mx 
~ ¢3a3(1 — e2) 
which is equal to the orbit-averaged precession frequency 
of the nodal axis. The expression given by Eq. (17) can 
be rewritten as (Peifker et al. 2020d), 


YL 3/2 = Our 5 (17) 


yxe [Ime 
2r, VV lt+e- 


Y= Qur =. (18) 
The Lense-Thirring precession rate of 54716 for the 
Ser A* spin of x = 0.5 is y, = 1.11 + 0.13 arcsec per 
year, which is larger than for $2 and $4711, see Peifsker 
et al. (2020d) for the comparison in terms of the FSS- 
population. 

Since $4716 appears to have one of the largest re- 
liable Lense-Thirring precession rates (S62 and $4714 


have large uncertainties), one can in principle turn the 
LT rate around to constrain the spin of Sgr A*. Ali 
et al. (2020) introduced the disk-like kinematic struc- 
ture of the S cluster in 3D space. Based on that, Fra- 
gione & Loeb (2020) constrain the spin of Sgr A* to 
x < 0.1 using the argument that the Lense-Thirring pre- 
cession would disrupt the disk structure (see also Fra- 
gione & Loeb 2022). $4716 seems to deviate from the 
disk structure with its inclination of ~ 161° being offset 
by more than 70° from the disk inclination peak close to 
90°. This could be caused by the vector resonant relax- 
ation (VRR), which is the fastest for the inner orbits of 
the S cluster. The timescale of the VRR behaves with 
the semi-major axis approximately as Typr x a7/? for 
the power-law slope y of the density cusp of late-type 
stars, n, x r~?. Therefore, the inclination distribution 
of S stars is expected to become wider towards smaller 
semi-major axes. $4716 has the maximum-likelihood 
longitude of the ascending node of 254716 ~ 150°. Us- 
ing this value, $4716 falls into the 2 distribution of the 
“black” disk with the mean value close to 180°. Under 
the assumption that the temporal evolution of 54716 is 
driven mainly by the Lense-Thirring precession, one can 
put the constraint on the Qur requiring that the Lense- 
Thirring timescale is comparable to or larger than the 
lifetime of S stars in order not to completely obliterate 
the non-isotropic 2) distribution of S stars. Using the 
expected deviation of 054716 from the mean value for 
the “black” disk, AQ = Q — OQga71¢ & 30°, we can write, 


AQ 
Tur = =— = Niife ~ 6.6 x 10®%yr. (19) 
Qur 
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The required precession rate is quite small, Qrr ~ 
0.0164 arcsec per year, which constrains the spin of 
Ser A* to y S 0.0074, i.e. consistent with y < 0.1 found 
by Fragione & Loeb (2020). The requirement of the well- 
ordered, non-isotropic 3D structure of the S cluster, i.e. 
consisting of multiple well-defined coherent disks in 3D, 
puts quite a stringent constraint on the Sgr A* spin to- 
wards small values. Even if the compact orbit of $4716 
has already made one turn of 2) by 7 radians during its 
lifetime, i.e., it would be now located where it was ini- 
tially within one of the disks, the precession rate would 
still be small Qur < 0.115 arcsec per year and the SMBH 
spin at most x < 0.052. 


4. DISCUSSION 


In this section, we categorize our findings and dis- 
cuss possible impacts on the detection of the compact 
and extended mass in the center of our galaxy. After 
the pericenter passage of $2, we reidentified $4711 in 
2020 with OSIRIS that were previously investigated in 
Peifsker et al. (2020d) and introduce a new S-star called 
S148. We would like to highlight a blend star event in 
2017 of S148 and $4716 that is discussed in detail in 
Appendix H. 


4.1. The observation of 54716 between 2003 and 2021 


While we continue to analyze the S-stars and the 
immediate environment, we recognized an unidentified 
stellar source in 2019 using public NIRC2 and NACO 
data. With our orbital tracking list, we identified the 
surrounding stars of this new stellar source neglecting 
an identification problem. Because of the proper mo- 
tion that is directed from east to west with a prominent 
curvature in 2019.30, 2019.36, 2019.49, and 2019.67, 
we were able to identify this new star in 2020.42 and 
2020.58 following a traceable trajectory. As GColl2022a 
claims to trace $29 by promoting a new orbital solution 
compared to Gillessen et al. (2009b), G17, GColl2021a, 
and Peifsker et al. (2021b), we provide an alternative in- 
terpretation of the interferometric GRAVITY data ob- 
served in 2021 with robust imaging data between 2003 
and 2020. 

By investigating published and archival data, we iden- 
tified this new source that we call $4716 in 16 epochs 
between 2003 and 2020 (see Table 8, Appendix A). Since 
crowding is a general problem in the S-cluster, 54716 is 
confused with $62 in 2004, 2011, and 2012. Moreover, 
the close distance of $2 hinders a confusion-free detec- 
tion of the star in 2016 (Appendix G). We underline 
the robustness of our finding by a photometric analy- 
sis covering data between 2003 and 2020. In agreement 
with the lo uncertainty range, we find an averaged K- 
band magnitude for $4716 of magx = 17.02+0.22 mag. 


Although fluctuations could be explained by a variable 
background, the extinction toward the Galactic center 
should also be considered (Fritz et al. 2011). Further- 
more, Schédel et al. (2010) reported a median extinction 
of Ax = 2.74+0.3 mag which is supported by the find- 
ings of Peifker et al. (2020c). Although it is plausible 
that the infrared magnitude of the general S-cluster pop- 
ulation can suffer from extinction variations, it is rather 
unlikely that these fluctuations show an imprint on the 
small $4716 orbit. The overall on-sky coverage of the 
S4716 orbit is in the order of around 0.6% compared 
to the S-cluster size. If the variation of the extinction 
causes the uncertainties observed for $4716, they are 
very well covered by our derived standard deviation of 
0.22 mag (for comparable values, see Schédel et al. 2010; 
Fritz et al. 2011; Gautam et al. 2019). We note that the 
faintest K-band magnitude values (Fig. 7) are observed 
with NIRC2/KECK implying a general characteristic of 
the detector. Although the wavelength coverage of the 
investigated data is comparable (see Table 8), the faint 
NIRC2 data points could imply a different detector sen- 
sitivity. Here, one has to take into account that the on- 
source integration time of the NACO data is increased 
compared to the KECK data (see Appendix A) result- 
ing in a presumably higher S/N ratio. Because faint 
sources like $4716 are sensitive to background emission, 
the different on-source integration time is a suitable ex- 
planation for the observed effect. We also note that we 
estimate the magnitude for $4716 with the ratio indi- 
cated in Eq. 10 which should be independent of the 
instrument. 

Furthermore, we witnessed a blend star event with a new 
identified S-star called $148 in 2017 (see Appendix H). 
In 2018, we reidentify $4716 at the expected projected 
on-sky position on its orbit. The data of 2016-2018 is 
shown in Appendix G where we provide photometric in- 
formation of $148. 

For the epochs 2007 and 2009, we extract the spectrum 
of $4716 using SINFONI data. The line-of-sight velocity 
of about -500 km/s and +1690 km/s matches the orbital 
fit model providing an independent confirmation of the 
detection. It is generally plausible that the chance of 
observing this star during the apoapse is increased due 
to the slow proper motion during this orbital section. 
In addition to the observation of $4716, we find $4711 at 
the expected position during its apoapse in 2020 using 
the orbital elements from Peifker et al. (2020d). Like- 
wise for 54716, the star $4711 can be easier observed 
during the apoapse. 


4.2. Systematics 
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Due to the number of different intruments, the anal- 

ysis is sensitive to systematic uncertainties. As we 
pointed out in Sec. 2.6, the data from different sources 
seem to result in agreeing high precission outcomes. 
However, Plewa et al. (2015); Plewa & Sari (2018) an- 
alyze the effects of the NACO detector on the data ob- 
served in GC. Plewa et al. furthermore discuss faint un- 
recognized objects that might influence the background 
emission and, therefore, the source determination. In, 
e.g., Peifker et al. (2020d), we faced these problems by 
identifying several new S-cluster members. Without a 
doubt, the S-cluster has many more sources as is known 
today (for a review, see Clénet et al. 2019). In addition 
to the analysis of the NACO detector and its influence on 
astrometric measurements, Jia et al. (2019) investigates 
positional uncertainties as a function of the magnitude 
of the K-band as observed with KECK. For a 17 mag 
faint star like $4716, the positional uncertainty is in the 
range of 0.1-1 mas. However, Jia et al. also show that 
uncertainties up to 10 mas are possible in the 17 mag 
K-band range. 
We can directly compare the position of $2/S0-2 with 
the orbital fit provided by Do et al. (2019) and 
GColl2022a. Furthermore, we choose $29 as a proper 
choice to investigate the dimensions of the systematics. 
This star was observed between 2002 and 2019 using 
mainly NACO data® (Peifker et al. 2021b). In Fig. 10, 
we find $29 in the expected position in two different 
epochs. We can conclude that the systematics and dif- 
ferent spatial resolutions are not the dominant origin of 
uncertainties. 


4.2.1. The proper motion of S29 


While we are not questioning the high precision mea- 
surements of GColl2022a, we would like to discuss a 
possible confusion in the presented data of the related 
publication. We remind the reader about the nomen- 
clature that is introduced in the introduction of this 
work: we distinguish between 529 (Peifker et al. 2021b), 
S29Gintessen (G17), and $29 (GColl2022a). While the 
first two listed publications provide different orbital so- 
lutions to the trajectory of S29 based on an almost 20- 
30 year data baseline, the common consensus is a peri- 
apse passage that does not match $29 as proposed by 
GColl2022a. Hence, we raise the question which star 
the authors actually have observed since the mentioned 
data baseline of almost 30 years analyzed in G17 covers 
a significant part of the S29@itlessen orbit. 
Unfortunately, the authors of GColl2022a missed the 
chance of indicating the precise epoch for the presented 


> Except for 2014 and 2019 where we used SINFONI data. 


S29 orbit. Since the S-stars are characterized by a high 
proper motion (Genzel et al. 2000), this will have an 
impact on the uncertainty of the here presented discus- 
sion. Anyhow, by analyzing Figure D.1. of GColl2022a, 
we extract x- and y- distances® of $29 with respect to 
Sgr A* which are listed in Table5. We use the values 


Year Axin [mas] Ay in [mas]  vVprop in [km/s] 


2015 148.1 237.0 - 
2016 140.7 222.2 562.1 
2017 125.9 177.0 1593.3 
2018 111.1 140.7 1353.9 
2019 88.9 88.9 1914.7 
2020 - - 1742.63 
2021 26.6 7.4 1742.63 


Table 5. Proper motion of S29 derived from GColl2022a 
(see their Fig. D.1.). The here listed proper motion is cal- 
culated by the year-to-year difference. The authors of the 
related publication do not provide a data point for 2020. 
Hence, we use the positional difference between 2019 and 
2021. It is implied, that the proper motion of $29 is at least 
in 2020 or 2021 too low. 


from Table 5 to derive the proper motion between two 
consecutive years by calculating the positional change 
over 365 days. We are aware that this is a rough ap- 
proximation, but it should serve the purpose since GC 
observations with NACO can not be executed through- 
out the whole year’. For the calculation of the associ- 
ated proper motion, we include the GC distance relation 
14mas = 100 AU = 150-10%km. With this, we calcu- 
late the related proper motion given in Table 5. For the 
data in 2020, the authors of GColl2022a are not provid- 
ing a data point. Therefore, we calculate the difference 
between 2019 and 2021 and divide the resulting value by 
2. For the 4 months in 2021 during the pericenter pas- 
sage of the star observed with GRAVITY, we calculate 
a proper motion of 6025 km/s. However, fitting a 3rd 
order polynomial to the data underlines the confusion 
regarding $29, S29, and $4716 (Fig. 11). As mentioned 
in GColl2022a, the star S29 is confused between 2015 
and 2017. We indicate in Fig. 11, that the disagreement 
with the orbital solution and the derived proper motion 
starts in 2017. We furthermore find that the star decel- 
erates in 2018, gets faster in 2019, and decelerates either 
in 2020 or 2021. We can conclude that the velocity of 


© We use the software Digitizer (Thomas Ott, MPE Garching) for 


the extraction of data points. 


7 From Paranal/Chile, the GC is only observable between end of 


February up until the beginning of October. 
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Figure 11. Proper motion of S29 as presented in 
GColl2022a. Because of the missing listening of the exact 
epoch per year in the related publication, we used a large 
x and y uncertainty of +0.5 years and +200 km/s. The 
red curve represents a 3rd order polynomial. Inspecting the 
proper motion underlines a identification problem for $29. 


the observed star in 2021.5 by GColl2022a exceeds the 
proper motion of $29 that is derived from the presented 
imaging data in the associated publication. 


4.3. Extrapolating high-precission data 


We would like to emphasize that we are not ques- 
tioning the capabilities of GRAVITY. However, using 
high-precision data and extrapolating these data to a 
crowded cluster may result in misidentifications. To un- 
derline this point, we shortly want to discuss the ob- 
servation of the linear motion of S62}incar (see Gillessen 
et al. 2009b, 2017). To date, it is not clear how a star 
can be gravitational bound to Sgr A* without showing 
any signs of acceleration or curvature. We discussed this 
already in Peifsker et al. (2021b) but we want to visual- 
ize the basic problem that emerges from this approach. 
For this, Fig. 12 shows the extracted data points from 
GColl2021a. Data prior to 2015 observed with NACO 
suggest a non-linear trajectory for S62)inear. Although 
we are aware of the noisy character of the NACO and 
SINFONI data, a non-linear trend for the data points is 
clearly observed before 2015. After 2015, the data ob- 
served with GRAVITY were extrapolated, resulting in 
the linear classification of S62. We point out that $55 
was classified as a linear star (Gillessen et al. 2009b) be- 
fore Meyer et al. (2012) derived a Keplerian orbit with 
an orbital period of 11.6 years. The misindetification of 
555 and $62 underlines our argument that extrapolating 
high-precision data into a crowded cluster may result in 
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Figure 12. Linear motion of S62 as proposed by Gillessen 
et al. (2009b), Gillessen et al. (2017), and GColl2022a. 
Data before 2015 was observed with NACO, after 2015 with 
GRAVITY. The non-linear motion of S62 is obvious for the 
data points before 2015 and follows the trend that we already 
presented in Peifsker et al. (2020a). The large gap between 
the NACO- and GRAVITY-data spans over 7 years. This 
data is extracted from GColl2021a. 


the confusion of stars. This confusion is reflected by the 
discontinuous velocity of $29 as shown in Fig. 11. 


4.4. Limitations of the analysis of the S-cluster using 
S300 


Inspecting the data sets reveals stars with a variable 
magnitude and high confusion. While these problems 
are already discussed in Sabha et al. (2012) and Peifsker 
et al. (2020a), we would like to discuss some of the chal- 
lenges with the robust observation of S62, S4711, and 
S4716 (see Fig. 10) over two decades. While all three 
listed stars can be observed in the majority of available 
data sets, S2 hinders the observation of $62, $4711, and 
54716 in the years during its periapse passage. As we 
know from the observation of S300 in GColl2022a, the 
diffuse background of the here presented imaging data 
can suffer from variations that are challenging to deter- 
mine (Sabha et al. 2012; Eckart et al. 2013). For $300, 
the range of K-band magnitudes are between 19.3-20.5 
mag. Based on the modeling of JaroszyNski (1999), it is 
plausible that other stars like $300 close to Sgr A* ex- 
ist that could interfere with the confusion-free detection 
of $62, $4711, and $4716. Donating such a star with 
S300*, we assume a K-band magnitude of 19.9 + 0.6 
mag. Using $2 as the reference star with a K-band mag- 
nitude of 14.15 mag, the (pixel-)counts from an obser- 
vation of $300* are 0.5% of the peak maximum. Since 
S$4711-S4715 are indicating the detection limit for the 
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here discussed imaging data, it is reasonable to subtract 
up to 2% of the image flux (i.e. background subtrac- 
tion). With this approach, the image information is pre- 
served since the flux of every observed source is above 
the upper limit of 2%. However, the 2% limit also shows 
that stars like $300* do exist inside the faint background 
and could presumably have an measurable impact on the 
extended mass. It is important to note that this 2% limit 
shows a different impact on the image area: the back- 
ground close to Sgr A* is higher compared to the outer 
region of the S-cluster (see Fig. 4 and Sabha et al. 2012). 
Hence, the chance of losing information especially close 
to the image edges is increased with the here discussed 
methods. 


4.5. The compact and extended mass 


From the MCMC statistics, we derive an enclosed 
Mass Mene Of (4.023 + 0.087) x 10°Me for the orbit 
of $4716 which numerically agrees with the results of 
Boehle et al. (2016). This derived result includes the 
compact Mass Mcom that is associated with Sgr A* and 
furthermore the extended mass m,,;. Here, Mmez+ iS a 
measure for the dark mass which should have a mea- 
surable imprint on the orbit of $4716. Deviations from 
the Keplerian elliptical orbit can arise either due to an 
extended mass inside a stellar orbit (retrograde peri- 
apse precession) or due to post-Newtonian effects of 
the curved spacetime around a point mass (prograde 
Schwarzschild precession, see Eq. (1); Weinberg 1972; 
Rubilar & Eckart 2001; Mouawad et al. 2005). In real- 
ity, both Newtonian and relativistic effects are expected 
to contribute to a smaller or greater extent. 

Following this approach, we define the measured en- 
closed mass as a sum of the compact and extended mass 
with 

Mence = Mext + Mcom - (20) 


When we use Eq. (20) in units of mene, we can write 
1 = Mext/Mence + Mcom/Menc Where as a consequence 
Meom = Menc yields Meee = O for a pure point-mass 
model, i.e. only the SMBH would contribute to the en- 
closed mass. For this, we assume a sphere with a volume 
Vsphere that includes both mass components. 

In the following, we assume that part of the mass located 
within the $4716 orbit is due to an extended mass Mezt 
component surrounding the otherwise heavy compact 
Mass Mcom Of SgrA*. The amount of an extended mass 
can be estimated by putting a limit on the correspond- 
ing retrograde Newtonian periastron shift, as long as 
it is large with respect to the prograde relativistic shift. 
The orbit of $4716 allows us to make a statement on the 
amount of mM. contained in the volume that is trans- 
versed by the $4716 star on its orbit. The distribution 


of any extended mass within Vspnere is unknown. How- 
ever, if Mex is more concentrated towards the SMBH or 
towards the rim of Vspnere, then the length of the orbital 
sections within that extended mass and correspondingly 
the amount of the Newtonian retrograde shift is dimin- 
ished. Hence, assuming the first-order homogeneous 
mass distribution puts a realistic upper limit on the 
amount of Meg:. Rubilar & Eckart (2001) as well as 
Mouawad et al. (2005) have investigated the effect of an 
extended mass on stellar orbits using the formalism pro- 
vided by Jiang & Lin (1985) in combination with orbit 
calculations using the post-Newtonian approximation of 
General Relativity (Weinberg 1972). For $4716, we will 
derive in the following a simple estimate of the order of 
magnitude of the extended mass within the Vspnere that 
is transversed along its orbit. 

For a body travelling on an elliptical orbit, the orbital 
time scale 7’, the semi-major axis a, and the central mass 
M are related via the third Kepler law 


3 GMT? 
a = 
Ar? 
Here, G is the gravitational constant. Furthermore, one 
finds for an elliptical orbit 


b=avl—-e? . (22) 


Here, 6 is the semi-minor axis and e the eccentricity 
of the orbit. If the amount of extended mass is Mezt, 
then varying the Mass from M to M + mez+ results 
in a change of orbital elements. Requiring that a for 
each semi-elliptical section of the rosetta-shaped New- 
tonian orbit remains constant, the changes are restricted 
to e and b, correspondingly. Hence, for successive semi- 
elliptical sections, the argument of periapsis w changes 
in a retrograde way. In the fortunate case of $4716, the 
orbit is almost face on, therefore the uncertainty of dw of 
the argument of periapsis w can immediately be trans- 
lated to an uncertainty of the compactness of the central 
mass M due to an additional extended mass component 
Mest: 
The MCMC fit of the $4716 orbit gives a lo uncertainty 
of dw = 0.023° (0.0004 rad). This corresponds to a pos- 
sible extended mass contribution of 7x 10° Me on the lo 
level and a 2x10* Mo contribution on the 3c level. This 
is more than 25 times smaller than the mass uncertainty 
given by the MCMC fit of the $4716 orbit. Our derived 
estimate of the extended mass is therefore on the same 
level as the robust estimate by Mouawad et al. (2005) 
for the 2.5 times larger $2 orbit and the corresponding 
limits summarized by Gillessen et al. (2009b) and G17. 
With mert < 2 x 104 Mo and mene ~ 4 X 10° 
Ma we can then directly derive a limit on the en- 
tire compact mass component associated with SgrA* as 


(21) 
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Meom < 3.98 x 10° Mo. By normalizing Eq. (20) to the 
entire enclosed mass Menc, we find Mezt/Mene S$ 5X10~3 
and Mcom/Mence < 0.995 while the sum of both terms 
equals unity. 


4.6. Could $4716 be a binary? 


Recent K-band observations with KECK suggest that 
the binary fraction in the Galactic center and, in par- 
ticular, the S-cluster is consistent with the local So- 
lar neighborhood fraction (Lefévre et al. 2009; Gautam 
et al. 2019). Moreover, Stephan et al. (2016) inferred a 
binary fraction of ~ 70% for the stellar population of the 
most recent star formation episode in the nuclear star 
cluster ~ 6 Myr ago (see e.g. Lu et al. 2013). Follow- 
ing the compelling evidence that stars form in multiple 
systems (Reipurth 2000; Looney et al. 2000; Sadavoy 
& Stahler 2017), we discuss the possibility that $4716 
could be a binary system. However, as we will show, 
the parameter space of such a system is limited in or- 
der to be stable enough. Because the largest variability 
amplitude would be achieved for an eclipsing system, 
we will focus our discussion on these systems. We note 
that the detectability of such a system in the S-cluster 
is at the limit of current instruments. This is underlined 
by the analysis of Gautam et al. (2019), who report an 
uncertainty of dmx, = 0.03 mag for the photometric 
study of the NSC. The authors found no evidence of an 
eclipsing binary in the S-cluster with a sample size of 17 
stars. Due to the limited sample size, the existence of 
an eclipsing binary system cannot be ruled out consider- 
ing an eclipsing binary fraction of approximately ~ 3% 
(Pfuhl et al. 2014; Gautam et al. 2019). 

For simplicity, we consider roughly equally luminous 
binary components. Then the variability amplitude in 
magnitudes is at most 6m < —2.5log(1/2) ~ 0.75 mag. 
Given the K-band magnitude uncertainty of 0.22 mag 
this is just at the limit of detection. 

For the pericenter distance of 54716 rp ~ 97.6 AU, the 
component semi-major axis can be at most, 


1/3 
Tp Mpbin 
AU 23 
97.6 oe & = » a 


dpin < 0.92 ( 


where the binary mass mpin is set to the inferred mass 
of $4716. For the component distance greater than apin, 
the binary would get disrupted at the periapse via the 
Hills mechanism. The corresponding orbital timescale 
of such a compact binary is 


2 dis 3/2 (my Ts 

Pa 180.9 (=) oe days. 

eine 0.92 AU (sae | a 
(24) 


If we consider the stellar diameter of D, = 2R, = 
2x 2.45 Re = 4.9Ro, the eclipse event would last at 


<S D,/vpin ~ 0.71 days, which is a dimming 
event too short to be practically detected with a sam- 
pling larger than one day. 

In other phases than the eclipse or for a non-eclipsing 
binary, components cannot be resolved out since the an- 
gular separation is < 0.11 mas, which is much smaller 
than the resolution limit of ~ 60 mas of eight-meter 
class telescopes. ‘The detection by spectral disentan- 
gling would also be challenging since the maximum ve- 
locity shift of spectral lines would be Umax ~ Upin ~& 
55.5kms~! for an edge-on system. It could potentially 
be larger for more compact systems though, however, 
the timescales of the velocity shift would also shorten 
(for the separation of 0.1 AU, the orbital timescale would 
be just 6.5 days) 

In summary, $4716 could only be a very compact bi- 
nary system with the separation of < 1AU. The de- 
tection of such a compact system is extremely challeng- 
ing in the crowded environment of the Galactic center 
since a potential eclipse event would last only ~ 0.7 days 
and the spectral line velocity shift would be at most 
~ 55.5kms~! for a given component separation. 


most tec < 


4.7. On the origin of the compact orbit 


The orbital compactness of $4716 with the semi-major 
axis of a = 1.94+0.02 mpc = 400 +4 AU is not straight- 
forward to explain in terms of in situ star-formation 
since the number density inside any cloud on these scales 
would have to be exceedingly large. This follows from 
the Roche density condition for the cloud to be tidally 
stable, 


M ‘a 
oche > 1.8 x 1025 : ( 
DE ee (sao 400 AU 


) ‘ cm? 
(25) 
Therefore the occurrence of a star on such a compact 
orbit requires a certain migration mechanism that brings 
the star from larger distances closer in. One possibility 
is the inward migration within a dense accretion disc 
due to “Type 1” torques (Goldreich & Tremaine 1980), 
which, however, works under the assumption that there 
has been a dense enough disc during the S-star lifetime. 
A simpler explanation for the inward migration of an 

S star is the occurrence of other bodies in the §S clus- 
ter that can effectively catalyze the semi-major axis and 
eccentricity changes. The Hills mechanism is the mecha- 
nism that involves a binary disruption. If the binary has 
the semi-major axis @pin and the mass Mpin, we obtain 
the tidal disruption radius that should coincide with the 
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54716 pericenter distance, 


Te =p = Opin 


Mpbin 
103 ( ai ) main)” AU. (26) 
1.4AU/ \ 10 Mo 


In case $4716 would be the captured star from the binary 
disruption, its mean expected semi-major axis would be 
(Hills 1991), 


mM, \2/3 
) ~4256AU. — (27) 


<A. >& 0.56apin ( 
Mpbin 


The corresponding eccentricity would be e, = 1—1"/ < 
ace >~ 0.976. Hence, these orbital elements differ from 
those of $4716 — the stellar orbit resulting from the bi- 
nary disruption is expected to be much more elongated. 
Even for the initial binary mass of mpin ~ 100 Mo and 
the semi-major axis of apin ~ 2.9 AU, the semi-major 
axis of the captured star would be < a, >~ 1899 AU 
and its eccentricity of e, = 0.95. 

A more likely way to shorten the semi-major axis of 
an 9 star is the three-body interaction within the S clus- 
ter — a close encounter of the star with a more massive 
perturber. If the original orbit is characterized by the 
orbital elements a, e, and 1, the close encounter with a 
more massive star or the intermediate mass black hole 
(IMBH) will change the orbital elements to a’, e’, and 
u’, while the Jacobi intergral or the so-called Tisserand 
parameter is preserved (Merritt 2013). If the perturbing 
body has the semi-major axis a,, the ‘Tisserand param- 
eter can be calculated as, 


1/2 
(4.64) = = +2 = - °)| COS L, (28) 


where vz is the inclination of an S star with respect 
to the perturber-SMBH plane. If the original orbit of 
54716 was nearly circular, then the close encounter has 
to occur for a ~ a, and a small inclination, hence 
T ~ 3. Considering the current orbital elements of 
54716, we can obtain a, from the condition T(a’, e’) ~ 3, 
where we assume that $4716 stays close to the orbital 
plane of the perturber. We obtain ap ~ 836 AU for 
a’ ~ 400 AU and e’ = 0.756. With the apocenter dis- 
tance of r, = a’(1+e’) ~ 702 AU, $4716 would approach 
the mean distance of the perturber as expected. Hence, 
the encounter of S stars with more massive bodies such 
as IMBHs could shorten their semi-major axes by at 
least a factor of two and make their orbits more eccen- 
tric. ‘The encounters could also be consecutive in an 
analogous way as comets interact with the large planets 


of the Solar system. This way $4716 could have mi- 
grated from the outer periphery of the S cluster all the 
way to its inner parts, though this needs to be verified 
in more detailed numerical simulations. 


5. CONCLUSION 


In this work, we present the observation of a new S- 
cluster member that we call $4716. This nomenclature 
follows the analysis presented Peifker et al. (2020d). We 
furthermore reidentify $4711 in 2020 using the OSIRIS 
science cam. For $4716, we derive an orbital period of 
about 4 years. Inspecting the 2020 imaging data under- 
lines the analysis of Peifker et al. (2020a), Peifsker et al. 
(2020d), Peifsker et al. (2021b), and this work, since we 
find within the uncertainties $62, $4711, $4716, and 
529 at the expected_position. We conclude that the 
periapse passage of $29 presented in GColl2022a is the 
result of the mixed observation of S62 and $4716. In 
the following, we list our key findings and results from 
the discussion section. 


1. $4716 can be observed with NACO, SINFONI, 
NIRC2, and OSIRIS between 2003 and 2020. The 
B8/9-V star orbits Sgr A* in 4 years at distances 
between 98 AU and 702 AU. The K-band magni- 
tude is about 17 mag and constant between 2003 
and 2020 within the lo uncertainty. 


2. $4716 is confused with $62 in 2004, 2011, and 
2012. In 2016, the observation of $4716 is ham- 
pered due to 82. Furthermore, we find a new 
S-star called $148 that is traceable over multiple 
epochs. In 2017, we observe a blend star event for 
$148 with $4716. This event is underlined by the 
photometric analysis for both stars between 2016 
and 2018. 


3. As a limit for the extended mass that may exist 
in addition to the compact mass of Sgr A* within 
the $4716 orbit, we derive 2x10* Mo (30). 


4. Calculated from the $4716 orbit, we find an en- 
closed mass of (4.023 + 0.087) x 10° Mo with a 
related distance of 8.028 kpc + 2% consistent with 
the literature (Akiyama et al. 2022). 


5. The data suggest that the claimed periapse pas- 
sage of the S-star $29 in 2021 by GColl2022a is 
instead the periapse passage of the here presented 
newly discovered S-cluster member $4716. 
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APPENDIX 


In this appendix, we will provide additional information about the analysis, follow up on the disputed identification 
of S29, and present supplementary data regarding the identification of S62 and $4716. We want to stress the observed 
blend star event in 2017. 

A. USED DATA 


Here, we list the used data for this work (see Table 6 & 7). For the KECK observations, we download the data from 
the KOA archive. The analyzed KECK data was observed at 2006-05-03 (21 exposures, 58.8 sec total integration time) 
, 2019-04-20 (14 exposures, 39.2 sec total integration time), and 2019-06-30 (10 exposures, 28.0 sec total integration 
time). We want to highlight that this data was also analyzed in Parsa et al. (2017), Peifker et al. (2019), Peifker et al. 
(2020a), Peifsker et al. (2020b), Peifsker et al. (2020c), Peifker et al. (2020d), Ali et al. (2020), Peifker et al. (2021a), 


Peifker et al. (2021b), and Peifsker et al. (2021c). 


NACO 
Date (UT) Observation ID number Total Seeing 
of exposure in arcsec 
exposures time(s) 
2003-06-13 071.B-0078(A) 102 510.00 0.6 
2005-07-25 075.B-0093(A) 330 343.76 1.0 
2005-07-27 075.B-0093(C) 158 291.09 0.9 
2005-07-29 075.B-0093(C) 101 151.74 0.9 
2005-07-30 075.B-0093(C) 187 254.07 0.9 
2005-07-30 075.B-0093(C) 266 468.50 0.9 
2005-08-02 075.B-0093(C) 80 155.77 0.9 
2007-03-04 078.B-0136(B) 48 39.86 0.85 
2007-03-20 078.B-0136(B) 96 76.19 0.85 
2007-04-04 179.B-0261(A) 63 49.87 0.95 
2007-05-15 079.B-0084(A) 116 181.88 0.75 
2008-02-23 179.B-0261(L) 12 86.11 0.9 
2008-03-13 179.B-0261(L) 96 71.49 0.8 
2008-04-08 179.B-0261(M) 96 71.98 0.8 
2010-03-29  183.B-0100(L) 96 74.13 0.9 
2010-05-09 183.B-0100(T) 12 16.63 0.85 
2010-05-09  183.B-0100(T) 24 42.13 0.85 
2010-06-12 183.B-0100(T) 24 A7.45 0.85 
2010-06-16 183.B-0100(S) 48 97.78 0.65 
2013-09-01  091.B-0183(B) 178 1780.00 1.05 
2016-03-22 594.B-0498(1) 144 6300 0.6 
2016-09-25 097.B-0216(B) 83 16.60 0.75 
2016-09-26 594.B-0498(K) 13 52.00 0.7 
2017-06-16 598.B-0043(L) 36 144.00 0.65 
2018-04-24 101.B-0052(B) 120 1200.00 0.55 
2019-09-27 5102.B-0086(H) Al 164.00 0.65 


Table 6. NACO K-band data observed between 2004-2019. Because several individual observations are stacked to a final 
mosaic, we list a range of the seeing for the corresponding used on-source exposures.Here, seeing is refering to the ESO archive 
keyword DIMM-Seeing-at-start. We note that observations with a seeing greater than 1.5 arcsec are excluded because the 
identification of individual sources is hampered. 
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Date Observation ID Amount of on source exposures Exp. Time Seeing 
Total Medium High 
(YYYY:MM:DD) (s) (arcsec) 
2013-04-06 091.B-0088(A) 8 0 8 600 Oar 
2013-04-07 091.B-0088(A) 3 0 3 600 0.6 
2013-04-08 091.B-0088(A) 9 0 6 600 1.3 
2013-04-09 091.B-0088(A) 8 1 t 600 0.7 
2013-04-10 091.B-0088(A) 3 0 3 600 0.9 
2013-08-28 091.B-0088(B) 10 1 6 600 0.7 
2013-08-29 091.B-0088(B) 7 2 4 600 0.9 
2013-08-30 091.B-0088(B) 4 2 0 600 0.9 
2013-08-31 091.B-0088(B) 6 0 4 600 0.7 
2013-09-23 091.B-0086(A) 6 0 0 600 1.3 
2013-09-25 091.B-0086(A) 2 1 0 600 Ld 
2013-09-26 091.B-0086(A) 3 1 1 600 0.7 
2014-02-27 092.B-0920(A) 4 1 3 600 0.8 
2014-02-28 091.B-0183(H) 6 3 1 400 0.9 
2014-03-01 091.B-0183(H) 11 2 4 400 0.8 
2014-03-02 091.B-0183(H) 3 0 0 400 0.9 
2014-03-11 092.B-0920(A) 11 2 9 400 1.0 
2014-03-12 092.B-0920(A) 13 8 5 400 0.9 
2014-03-26 092.B-0009(C) 9 3 5 400 0.9 
2014-03-27 092.B-0009(C) 18 7 5 400 1.0 
2014-04-02 093.B-0932(A) 18 6 1 400 0.9 
2014-04-03 093.B-0932(A) 18 1 17 400 0.8 
2014-04-04 093.B-0932(B) 21 1 20 400 0.8 
2014-04-06 093.B-0092(A) 5 2 3 400 0.7 
2014-04-08 093.B-0218(A) 5 1 0 600 1.1 
2014-04-09 093.B-0218(A) 6 0 6 600 0.8 
2014-04-10 093.B-0218(A) 14 4 10 600 0.8 
2014-05-08 093.B-0217(F) 14 0 14 600 0.7 
2014-05-09 093.B-0218(D) 18 3 13 600 0.8 
2014-06-09 093.B-0092(E) 14 3 0 400 0.8 
2014-06-10 092.B-0398(A) /093.B-0092(E) 5 4 0 400/600 0.7 
2014-07-08 092.B-0398(A) 6 1 3 600 0.7 
2014-07-13 092.B-0398(A) 4 0 2 600 nan 
2014-07-18 092.B-0398(A) /093.B-0218(D) 1 0 0 600 1.1 
2014-08-18 093.B-0218(D) 2 0 1 600 1.1 
2014-08-26 093.B-0092(G) 4 3 0 400 1.3 
2014-08-31 093.B-0218(B) 6 3 1 600 0.9 
2014-09-07 093.B-0092(F) 2 0 0 400 0.8 
2015-04-12 095.B-0036(A) 18 2 0 400 0.9 
2015-04-13 095.B-0036(A) 13 7 0 400 0.9 
2015-04-14 095.B-0036(A) 5 1 0 400 0.9 
2015-04-15 095.B-0036(A) 23 13 10 400 0.9 
2015-08-01 095.B-0036(C) 23 i 8 400 1.3 
2015-09-05 095.B-0036(D) 17 11 4 400 1.5 


Table 7. SINFONI data observed in 2013, 2014, and 2015. From the H+K data cubes, we extract the continuum images 
covering the wavelength range of 2.0 — 2.2 4m. We furthermore indicate the related seeing for every observation night. Because 
seeing changes for every single observation, we provide an averaged value that represents the overall quality of the observation 
night. We note that also under bad conditions (seeing larger than 1.5 arcsec), high-quality observations are, in principle, possible. 
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In Table 8, we indicate the instruments related to the data points shown in Fig. 6 and Fig. 9. Although the filters 
used in this analysis cover slightly different spectral ranges, the spectroscopic offset is rather small in terms of absolute 
numbers. Because we assume that the magnitude ratio between stars is constant for the investigated sample, the 
absolute sensitivity should not be a major contributing factor to the uncertainties. Since we use $2 as a reference 
star with a related magnitude of 14.1 mag in the K-band, the ratio of another star should be independent of the used 
instrument if the frame (sky correction, on-source integration time, detector health) of the observations is comparable. 
Consult Sec. 3 for an additional discussion. 


Epoch AR.A. (as) MADEC. (as) AR.A. error (as) MADEC. error (as) A [um] Resolution (as) Instrument 


2003.44 -0.041 0.066 0.007 0.007 1.97-2.32 0.039 NACO 
2005.26 -0.033 -0.002 0.007 0.007 1.97-2.32 0.046 NACO 
2005.58 0.019 0.033 0.007 0.007 1.97-2.32 0.039 NACO 
2006.33 -0.015 0.077 0.007 0.007 1.94-2.29 0.029 NIRC2 
2007.25 -0.039 0.073 0.007 0.007 1.97-2.32 0.030 NACO 
2008.46 -0.053 0.045 0.007 0.007 1.97-2.32 0.033 NACO 
2010.24 -0.017 0.079 0.013 0.013 1.97-2.32 0.039 NACO 
2013.40 0.006 -0.013 0.007 0.007 2.00-2.20 0.037 SINFONI 
2013.66 0.013 0.000 0.007 0.007 1.97-2.32 0.046 NACO 
2014.50 -0.024 0.070 0.007 0.007 2.00-2.20 0.037 SINFONI 
2015.30 -0.050 0.062 0.007 0.007 2.00-2.20 0.042 SINFONI 
2019.30 -0.027 0.073 0.007 0.007 1.97-2.32 0.029 NIRC2 
2019.67 -0.044 0.064 0.005 0.005 1.94-2.29 0.028 NIRC2 
2019.73 -0.053 0.066 0.007 0.007 1.97-2.32 0.038 NACO 
2020.42 -0.050 0.050 0.005 0.005 2.12-2.22 0.039 OSIRIS 
2020.58 -0.050 0.050 0.005 0.005 2.12-2.22 0.029 OSIRIS 
2021.24 -0.051 0.034 0.0005 0.0005 2.00-2.40 65 x 10°°* GRAVITY 
2021.41 0.0029 -0.01176 0.0005 0.0005 2.00-2.40 65 x 10°°* GRAVITY 
2021.47 0.0129 -0.0029 0.0005 0.0005 2.00-2.40 65 x 10°°* GRAVITY 
2021.56 0.0182 0.00941 0.0005 0.0005 2.00-2.40 65 x 10°°* GRAVITY 


Table 8. List of derived positions and the corresponding uncertainty for $4716. Using a GAUSSIAN fit to derive the position 
of $4716 is accompanied by an underestimated uncertainty that does not reflect the noise character of the region (Peifker et al. 
2021c). An uncertainty of +5 — 7 x 10~%as reflects positional uncertainties in combination with crowding challenges in the 
S-cluster. In addition, we indicate the resolution of the related data by measuring the FWHM of 82. For GRAVITY data, we 
use the astrometric accuracy indicated by GColl2022a. 


B. MCMC STATISTICS 


Here, we present the outcome of the MCMC simulations (Fig. 13). A compact distribution indicates a well-defined 
input parameter. For transparency, we are aware of the second peak for the inclination close to i = 2.90 = 166.1°. 
Of course, this implies that other solutions to the Keplerian fit simulations may also produce satisfying results. 
Nevertheless, we would like to guide the reader to Fig. 9 where we show the Keplerian solution for the observed 
data points for $4716. Since the Keplerian solution reasonably matches the observations, we use the result for the 
inclination from the MCMC statistics to determine a proper uncertainty range. For the uncertainty, we use the mean 
error WIth Imsawen = 166.1°—161.2° = 2.4°. With a mean value for the inclination of imean = 163.6°, we estimated a 


final value of 163.6° + 2.4°. 
C. MASS AND DISTANCE OF SGR A* 


Since we executed the MCMC simulations with iteration steps of i=10000 several times, we noticed a strong correla- 
tion of the mass and distance for Sgr A*. The results of these runs do show agreeing values indicating that $4716 is a 
perfect candidate for the evaluation of the enclosed mass and the distance of the SMBH. In the following, we present 
the related figures of the MCMC runs (Fig. 14-16). 
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Figure 13. MCMC simulation of the orbital elements given in Table 2 for $4716. Overall, we find a convincing agreement with 
the input parameters which is reflected in the compact shape of the shown output. 


D. LINE OF SIGHT VELOCITY OF $4716 


Here we show the related LOS velocity to the spectrum shown in Fig. 8. The fit presented in Fig. 17 is purely based 
on the outcome of the Keplerian orbit presented in Fig. 6. This means, we derive the orbital solution and the related 
LOS evolution along with the orbit independently from the estimated Dopplershifted velocity of almost -500 km/s in 
2007 and 1690 km/s in 2009 (Fig. 18). With this approach, we have an independent parameter to verify the orbital 
solution and the related source detection for these epochs. 
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Figure 14. Results of the first MCMC run. As described in the text, we use 10* iteration steps. Because of the open boundaries, 
the range of possible solutions for the mass and the distance covers almost 10% of the peak value. 
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Figure 15. Results of the second MCMC run. Like the results presented in Fig. 14, we use 10° iterations steps for the MCMC 
simulations. However, here the uncertainty range is decreased to 5%. 
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Figure 16. Results of the third MCMC run. Here, the uncertainty range is about 10% like the results presented in Fig. 14. 
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Figure 17. Line of sight velocity for 54716. Here we used SINFONI data of 2007 and 2009 to derive the Doppler-shifted LOS 
velocity of $4716 with the Bry absorption line. The green data points are not fitted and represent an independent parameter 
to confirm the here presented analysis. 


E. COMMENT ON THE DETECTION OF S62 IN 2019 


While we highly appreciate the attempt of observing $62 with GRAVITY in 2019 GColl2021a, we would like to 
address some challenges regarding the interpretation of the presented data. In Tab. 9, we present the projected 
positions of S29@ittessen (G17), S62 (Peifsker et al. 2020a), and the GRAVITY star that we denote with $29 observed 
in 2019 and analyzed in GColl2021a. Considering the close distance of S62 to $29 in 2019 results in a challenging 
interpretation, to connect the latter star to S29Gitlessen. To visualize this point, we created a mock image with the 
FOV of GRAVITY and the positions of the three stars in 2019.66 and 2019.75 as listed in Table 9. The authors of 
GColl2021a argue that the proper motion of the observed GRAVITY star does not match $62 but rather S29Gittessen- 
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Figure 18. Spectrum of $4716 in 2009 as observed with SINFONI. In contrast to 2007 (Fig. 8), the velocity is redshifted with 
about 1700 km/s with respect to the Bry restwavelength at 2.1661um. 


The authors ignore the predicted periapse passage for S29Giljessen but more importantly, the increasing distance to the 
observed GRAVITY star. Compared to $62, S29Gitlessen is over 400% away from the GRAVITY star. Therefore, we 
propose that S62 = S29 = GRAVITY star in 2019. 


F. THE DETECTION OF S62 


In this section, we show the same data as presented in Fig. 5 but with the $62 orbit. Even though we have presented 
numerous observations of this star (Peifker et al. 2020a,d; Peifsker et al. 2021c), we revisited the available data with 
the approach of using a combined PSF (APSF -+airy ring model) for the Lucy Richardson deconvolution algorithm. As 
expected, the data shown in Fig. 20 clearly shows the star S62 at the expected position on its orbit around Ser A*. 


G. OBSERVATION OF $4716 IN 2016, 2017, AND 2018 


With a K-band magnitude of 14.15 mag, the bright S-cluster member $2 is the dominant source in the cluster 
(Schédel et al. 2002; Habibi et al. 2017). We already faced the observational influence of S2 regarding S62 in 2015 
(Peifsker et al. 2020a). For the same epoch, the projected position of S2 comes close to the apoapse of $4716 (Fig. 
5 and Fig. 20). Therefore, we will investigate the gap between 2016 and 2018 for the possible detection of 54716 
besides the influence of S2 (see Fig. 21). We would like to note that the S-cluster is constantly monitored. Due to 
blending, confusion, variable background, and general data quality related challenges, the number of reported stars 
in the literature (see, e.g., G17) is much lower than the observed stars (see Peifker et al. 2020d). Therefore, we will 
introduce a new S-star that crosses the orbit of $4716. We call this new identified star $148, following the widely 
accepted nomenclature (Ali et al. 2020). In Sec. H, we will elaborate on the analysis of S148. We find that $4716 is 
confused with $148 in 2017 (see the circled star in Fig. 21 in 2017.4). The K-band magnitude of 5148 is almost 1 
mag higher in 2017 (16.0 mag) compared to 2016 (16.9 mag), underlining the confusion with $4716. In contrast, the 
circled star in 2018.3 (Fig. 21) shows a magnitude of 17.0 mag, indicating the reidentification of $4716 (see Sec. 3.2) 
in this epoch. This followed by the presented identification in 2019 (Fig. 5) and 2020 (Fig. 10). 


H. THE NEWLY IDENTIFIED S-CLUSTER STAR $148 


Here, we will present the orbital elements of $148. In 2016, S148 passed Sgr A* (see Fig. 21) and is located close 
(about 12.5-25 mas) to S63. In agreement with the presented list for verifying newly identified S-stars (Sec. J) we 
find this star in 3 consecutive years on a Keplerian orbit at the expected position in 2016.5, 2016.7, 2017.4, and 2018.3 
using the LR algorithm. In 2018.3, we additionally observe $148 using a SM-filter (see Sec. 2). 

The orbital elements are displayed in Tab 10, the solution for the Kepler fit is shown in Fig. 22. The uncertainties 
are adapted to account for the limited orbit coverage of the data because of confusion with $2 before 2016 and after 
2018. As mentioned in the previous section (Sec. G), the star shows a variable magnitude between 2016 and 2018. 
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Figure 19. Visualization of the source confusion as summarized in Table 9 in 2019.66 and 2019.75. Here, the green dashed 
circle indicates the fiber FOV of GRAVITY where the authors of GRAVITY Collaboration et al. (2021) claim the observation 
of a stellar source that is located in the center of these images. As it is implied by this figure, the closest source to the center 
of the GRAVITY fiber FOV is S62. The predicted position of S29Ginessen based on compromised orbital elements taken from 
G17 (see clarifying comments in Peifsker et al. 2020a) is at too large a distance from the center of the GRAVITY fiber FOV in 
order to still be considered. It is important to note that no bright source was located at the predicted S29Gitessen in 2019.66 
and 2019.75. Hence, we use a lime and black cross hatched symbol. Consequently we are allowed to raise the question, why 
the authors of GRAVITY Collaboration et al. (2021) came to the conclusion of neglecting the confusion free identification of 
S62. This misidentification is a key element that gave rise to the $29 orbit as published now in Gravity Collaboration et al. 
(2022a). We would like to express our doubts about this identification since the authors pointed GRAVITY to the position 
of S62 based on our proposal (see Peissker et al. 2019) and private communication to the GRAVITY observing team on the 
Paranal mountain. Here, the figure size is 150x150 mas, the diameter of the dashed green circle is about 74 mas. 


For the data of 2016, we measure a K-band magnitude of 16.9+0.1 mag for $148. The star gets significantly brighter 
by almost 1 mag in 2017, resulting in a K-band magnitude of 16.0+0.1 mag. This is followed by a decreased K-band 
magnitude to 17.0+0.1 in 2018 for $148. We can safely conclude that $4716 combined with $148 creates a so-called 
blend star (Sabha et al. 2012; Eckart et al. 2013). This apparent star is the result of at least 2 fainter stars and shows 
a higher K-band magnitude compared to the individual S-stars $148 and $4716. The created blend star is an example 
for the checklist introduced in the next section and underlines the need for longer data baselines. A detailed analysis 
of this event and a comprehensive analysis of this star exceeds the scope of this work. However, we will study $148 
and the corresponding blend star event in detail in an upcoming publication (Peifker et al., in prep.). 


I. MASS RELATION IN THE S-CLUSTER 


Here we outline an accessible method for deriving the mass of the stellar members of the S-cluster. For this, we 
use the estimated mass and K-band magnitudes of Habibi et al. (2017). The authors of Habibi et al. use a 12 year 
data baseline achieved with SINFONI. We follow the inspired analysis of Peifsker et al. (2020d) and fit a 1 dimensional 
polynom to the derived mass and magnitude of the S-stars (Fig. 23). We use log ars = m-x+b where x denotes the 


K-band magnitude mx. With the data from Habibi et al. (2017) and Fig. 23, we find log its = 0.1925 - mx + 3.885. 
J. CHECKLIST FOR THE DATA 


Here we will list some key aspects to verify findings with the here presented approach using high-pass filter for the 
analysis. Of course, this list can be applied to previous and future publications. 
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Projected positions of S29Gitlessen, 562, and $29 in 2019 

Epoch S29cinessen x S29cittesseny S622 S62, S29, $29, 
2019.66 -98.6 +108.6 -85.3 +97.2 -87.6 +92.5 
2019.75 -97.3 +105.3 -84.4 +96.7 -85.9 +88.7 


Table 9. Positions of S29Gintessen, S62, and $29 observed in 2019. We extract these values from the related publications 
(GColl2021a, GColl2022a) and show a visualization of this table in Fig. 19. We would like to note that there is a striking 
distinction between the positions of S29Gitlessen and $29. 
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2015.3 2019.3 2019.5 2019.7 


Figure 20. Detection of S62 in the here presented data set. The observations are carried out with the VLT (NACO and 
SINFONI) and the KECK (NIRC2) telescope. In 2015, S62 is confused with S2. Likewise in Fig. 5, the data observed with 
NIRC2/KECK for the epochs 2006.3, 2019.3, and 2019.5. We furthmore use SINFONI/VLTA data for the epochs 2013.4, 2014.5, 
and 2015.3. The remaining epochs are observed with NACO/VLT. Magnitude variations do not the represent the intrisic flux 
of the star. These variations are related to the background, data quality, nearby sources, and the contrast (which is choosen to 


show the surrounding single stars and $62). A detailed analysis of the K-band magnitude of S62 can be found in Peifsker et al. 
(2021b). 
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2016.7 2017.4 2018.3 


Figure 21. K-band view on the environment of Sgr A* (indicated by a x) observed with NACO showing $2 and $4716 for 
2016.7, 2017.4, and 2018.3. In 2016, we illustrate the expected position of $4716 in the related epoch with a pink colored patch 
on the green orbit. The $2 orbit is shown in lime. The yellow trajectory shows the orbit of $148. North is up, east is to the left. 


Source a (mpc) e€ i(°) w(°) Q(°) taicsent (VP) 
S148 31.13 + 0.1 0.971 + 0.05 111.154 5.00 174.75 + 5.00 359.24 + 5.00 2016.8 + 0.2 


Table 10. Orbital elements of $148. The uncertainties are an upper limit and reflect the poor data baseline showing the 
position of $148. A complete analysis of $148 would exceed the scope of this work. However, we find that $148 is following the 
disk structure of the S-cluster as proposed by Ali et al. (2020). 


1. Independent observations: if possible, one should use several instruments /telescopes. 


2. Long term: to avoid blend star scenarios, at least 3 consecutive years of observations should exclude this possibility 
because of the high proper motion of the S-cluster members. 


3. Spectrum: the Kepler fit arranges a possible orbit in the 3 dimensional space based on data points observed on 
the projected sky. Hence, the result of the fit also delivers a solution for the LOS velocity. If possible, a spectrum 
of the observed star can be used to verify the finding as an independent parameter. 


4. Inside a reasonable uncertainty range, the proper motion should coincide with the observation of the star. 
5. The latter bullet point can also be applied to the magnitude. 


6. Because of the stellar density of the cluster, unidentified stars can be located at the airy rings of known S-stars. 
If the star can only be detected on the airy rings, it is reasonable to assume that the finding is an artefact. 
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Figure 22. Orbital solution for $148 based on the measurements in 2016.5, 2016.7, 2017.4, and 2018.3. Because of the increasing 
distance between $148 and 82, this solution is likely to be improved in the near future. 
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Figure 23. Relation between log(M/Mo) and the K-band magnitude mx. We use the estimated mass and magnitude of Habibi 
et al. (2017). 


